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The critical behaviour of 3-dimensionaI disordered systems in a strong magnetic field is investi- 
gated by analysing the spectral fluctuations of the energy spectrum. The level spacing distribution 
P(s) as well the Dyson-Mehta statistics A3(L) are considered. Except for the small s behaviour 
of P{s), which depends on the presence or absence of a magnetic field but not on its strength, the 
other quantities, large s behaviour of P{s) and A3(L), turn out to be independent of its presence 
or absence and therefore of the time reversal symmetry. This suggests that the metal-insulator 
transition, which is defined by the symmetry of the system at the critical point, is independent of 
the presence or absence of the time reversal symmetry. 
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It is now well know that 3-dimensional disordered sys- 
tems exhibit a metal-insulator transition (MIT) as a func- 
tion of the disorder. Although a lot of work has already 
been devoted to this problem |Q] many points still remain 
unclear. One of them is the influence of the magnetic field 
at the MIT. In 2-dimensional systems the magnetic field 
has a dramatic effect. Whereas no MIT is expected with- 
out a magnetic field it turns out that systems in a mag- 
netic field exhibit one. Delocalised states appear at the 
centre of the Landau bands and quantities like the critical 
disorder or the critical exponent can be computed ||^ . In 
3-dimensional systems it was assumed, although there is 
no such dramatic effect, that the critical behaviour at the 
MIT can be cast into three different universality classes 
according to some general symmetries of the system: or- 
thogonal (with time reversal symmetry), unitary (with- 
out time reversal symmetry, e.g. with a magnetic field) 
and symplectic (with spin orbit-couphng) [Q. One then 
expects different critical exponents related to the MIT 
for the three different universality classes. Surprisingly, 
in spite of the change of universality class, the same value 
of the critical exponent has been found, numerically, with 
and without a magnetic field j4|j|] . A convenient way to 
study this problem is to resort to random matrix the- 
ory (RMT) and energy level statistics (ELS). It has been 
shown in the case without magnetic field, that be- 
side the two expected statistics, namely the Gaussian or- 
thogonal ensemble (GOE) for the metallic regime and the 
Poisson ensemble (PE) for the insulating regime, a third 
statistic, called the critical ensemble (CE), occurs only 
exactly at the critical point. It then seems natural that 
the universality class of the system should be determined 
by the symmetry of the system at the critical point with 
the results quoted above being interpreted as the 
sign of the independence of the universality class, at the 
critical point, to teh breaking of time reversal symmetry. 



This has been checked in the case of a system subject 
to an Aharonov-Bohm (AB)-flux Despite the fact 

there is no magnetic field inside the system the AB-fiux 
breaks the time reversal symmetry and one observes a 
transition this time from the Gaussian unitary ensemble 
(GUE) to the PE. But the CE, the critical exponents, 
does not depend on the time reversal symmetry and so 

A reservation which has to be considered with this kind 
of system is that a change of the vector potential is, by 
means of a gauge transformation, equivalent to a change 
of the boundary conditions and therefore should have no 
influence in the localised regime. Despite this it must 
be stressed that, when using periodic boundary condi- 
tions, the vector potential cannot be gauged away and 
the system is not equivalent to a system invariant under 
orthogonal transformation. Nevertheless, in this paper 
we would like to address the problem now with a mag- 
netic field inside the system and to show that the conclu- 
sions derived from the case with the AB-flux still remain 
valid. 

In order to investigate the MIT with a magnetic field 
we consider the usual Anderson Hamiltonian, 



(1) 



where the sites i=lmn are distributed regularly in 3D 
space, e.g. on a simple cubic lattice. Only interactions 
with the nearest neighbours are considered. The site en- 
ergy eimn is described by a stochastic variable. In the 
present investigation we use a box distribution with vari- 
ance (T^ = W^/12. W represents the disorder and is the 
critical parameter. Introducing a homogeneous magnetic 
field B = i?( 1,1,1) leads to complex hopping elements 
(Peierls factors) Using the following vector poten- 
tial A = B{z, X, y) and imposing \Vimn,i'm'n' | = 1, which 
defines the energy scale, yields for Vimn,i'm'n'- 
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exp(±i^2i2^) if /' = / ± 1, to' = TO, n' = n 
exp(±ii ^^°^^" ) if r = ;,TO' = m±l,n' = ?i (2) 
exp(±z ^^° if r = ;,TO' = TO,n' = n±l , 

where a is the step of the lattice, a — is the the 

number of flux quanta h/e per plaquette. The boundary 
conditions are taken to be periodic. 

Based on this Hamiltonian, the MIT in the presence 
of a magnetic field will be studied by the ELS method, 
i.e. via the fluctuations of the energy spectrum Be- 
fore giving the results we shortly review the ELS method. 
Starting from Eq.(l) the energy spectrum was computed 
by means of the Lanczos algorithm (which is suited to di- 
agonalise such very sparse secular matrices) for systems 
of size M X M X M with M = 13, 17, and 21, disor- 
der W ranging from 3 to 40 and phase a from 0.005 to 
0.250. The number of different realizations of the random 
site energies eimn was chosen so that about 10^ eigenval- 
ues were obtained for every pair of parameters (M, W) 
which means between 25 and 90 realizations, for which 
only half of the spectrum around the band centre is con- 
sidered so that the results do not deteriorate due to the 
strongly localised states near the band edges. After un- 
folding the spectrum obtained, the fluctuations can be 
appropriately characterised [ p| by means of the spacing 
distribution P(s) and the Dyson-Mehta statistics A3(L). 
P{s) measures the level repulsion, it is normalised, as 
is its first moment, because the spectrum is unfolded. 
A^^L) measures the spectral rigidity. 



case with time reversal symmetry , there is no longer 
one point where all the curves intersect but rather one 
which seems to appear only for the curves in the metallic 
regime. 

The next step is now to find where the MIT takes place. 
For this one uses the fact that the quantities we are con- 
sidering here, P{s) and A3(_L), depend on the system size 
except at the critical point where they are scale invari- 
ant. This is directly related to the fact that the MIT is 
a second order transition and that finite-size scaling laws 
will apply close to the transition. From a numerical point 
view it has been shown that A3(L) gives more accurate 
results. So we calculate S{M,W) = 5^/0 A.3{L)dL as 
a function of M and W, the critical disorder, Wc, being 
given by W for which S{M, W) is independent of M [ p^ . 
The results are shown in Fig. 0. 





FIG. 1. Spacing distribution, P(s), for M = 21 and dif- 
ferent disorders W . The full curves give the asymptotic 
regimes namely, GUE for the diffusive regime and PE for 
the localised regime. We see the transition between GUE 
and PE as a function of the disorder W . 

In Fig. the results for the spacing distribution P{s) 
are reported. One finds, as expected, the GUE and the 
PE regimes for small and large disorder respectively as 
well as the transition between them as a function of W. 
But less expected is the fact that, in contrast to the 



FIG. 2. 5 as a function of M and W". o : M = 13, □ : 
M — 17 and o : A/ = 21. The critical disorder is given by 
the point for which S is independent of M. Wc ranges from 
16.5 without a magnetic field up to 18.75 where Wc seems 
to become independent of the strength of the magnetic 
field. 

In contrast to the results obtained with an AB-flux 
Ip^ , Wc depends on the vector potential. This is due to 
the fact that now we have a magnetic fleld inside the sys- 
tem which breaks the weak localisation phenomena [0. 
It is then more difficult to localise the states which means 
an upwards shift of Wc- Wc ranges now from 16.5 with- 
out magnetic field up to 18.75 where Wc seems to become 
independent of the strength of the magnetic field. A sim- 
ilar behaviour has already been observed by Henneke et 
al g. With Wc it is now possible to study P{s) and 
A3(i) at the critical point. The first feature we would 
like to consider is the small s behaviour of P{s). In fact 
for numerical reasons what we calculate is the cumula- 
tive level-spacing distribution I{s) = P{s')ds' which 
allows a better study of the small s behaviour of P{s). 
The results are plotted in Fig. |^. As we can see the 
curves depend on the absence or presence of a magnetic 
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field but are independent of its strength. 
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FIG. 3. In-ln plot of the cumulative level-spacing dis- 
tribution 7(s) for small s. The curves depends on the 
absence or presence of a magnetic field but are indepen- 
dent of its strength. In both cases the s behaviour is the 
same as in the metallic regime but the constants in front 
of s and are higher, showing a decrease of the level 
repulsion. 

One finds, for small s, P{s) oc as^ with magnetic field 
and P(s) oc hs without. The behaviour is the same as 
in the metallic regime except for the pre-factors a and 
6 which are now both higher showing a decrease of the 
level repulsion. This comes from the nature of the criti- 
cal states which are multi- fractal . The overlap of the 
wave functions between the different sites, although still 
present, is now much weaker than in the metallic regime 
which decreases the level repulsion between the states. 




FIG. 4. Large s behaviour of P(s). The curves are no 
longer only independent of the strength of the magnetic 
field as for small s but are independent of the presence or 
absence of the magnetic field. In the inset exp —As, with 
A ~ 1.8, (Fit 2) and exp-(A's Cs^'^''), with A! ~ 1.37 
and C ~ 0.26, (Fit 1) have been fitted to the tail of P(s). 



In Fig. |j the other limit, namely the large s behaviour 
of -P(s), is reported. The curves are no longer only inde- 
pendent of the strength of the magnetic field, as for small 
s, but are even independent of the presence or absence of 
the magnetic field. The question of the shape, for large 
s, of P(s) is rather controversial. It was first claimed |^ 
that P(s) should be oc exp — As. Then following some 
new analytical calculations at the critical point Aronov 
et al. ^ proposed that P{s) = Bs'^ exp(-A/3s2-T) with 
7 = 1 — 1/vd where v is the critical exponent of the 
correlation length £^oo- It has to be noted that some 
doubts have been recently raised concerning the valid- 
ity of these results jl^. Numerically the situation is 
not really clearer. When one tries to fit a curve like 
P{s) ~ Bs exp(— As'') for the case without magnetic field 
one finds an 77 between 1.20 and 1.30 in a good 

agreement with the results obtains by Aronov et al. . 
But when one looks just at the tail of P{s) one can see 
some discrepancies between the fit and the numerical re- 
sults. Following some other numerical calculations it was 
then claimed ||l^ that P{s) was indeed oc exp — As as first 
suggested But the fit was done only for the large s 
and not on the whole curve as we did previously [p|,p^. 
In the inset of Fig. ^ we show that, when considering 
only the tail of -P(s), exp — As with A ~ 1.8 (Fit 2) gives 
indeed a good result but so does exp -{A'sA-Cs^-"^^) (Fit 
1). It is therefore very difficult to draw a definitive con- 
clusion about the problem just based on the fit of the tail 
of P{s). In the case with a magnetic field it turned out 
to be very difficult to fit a whole curve using the result 
of Aronov et al. which could be related, somehow, to 
the absence of a fixed point in P{s) as noticed above. It 
clearly shows that more work needs to be done in order 
to understand properly the shape of P{s) at the critical 
point. 
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FIG. 5. Dyson- Mehta statistics (A3(L)) at the critical 
point for various magnetic fields. As for the large s be- 
haviour of P(s), ^z[L) is not only independent of the 
strength of the magnetic field but of its presence or ab- 
sence too. 
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Finally we compare A3(L) for the different magnetic 
fields at the critical point in Fig. |[ As for the large s 
behaviour of P{s), A3(L) is not only independent of the 
strength of the magnetic field but seems to be indepen- 
dent of its presence or absence too. The shape of ^^{L) 
contains a term linear in L as well as a non-linear term 
L'^ , with a; < 1 and has already been studied in a previ- 
ous work In the inset of Fig. ^ we show how A^^L) 
varies, when we slightly change W, to give an idea about 
the accuracy of the curves keeping in mind that the un- 
certainty for Wc is ±0.25. These results suggest that, 
although the magnetic field brings some modifications at 
the critical point, as for the critical disorder Wc or the 
small s behaviour of P(s), the symmetry of the system 
is indeed independent of the magnetic field and then of 
the time reversal symmetry. What happens is, when the 
system reaches the critical point, for increasing disorder, 
the 0{N) or U(iV) symmetry is spontaneously broken. 
The reason for this break is obvious when one considers 
the localised regime. If the matrix would be O(A^) or 
U(A^) invariant it would then be possible to construct 
delocalised states by linear combinations of the localised 
ones which is clearly not allowed. In RMT this phenom- 
ena has led to the concept of a preferential basis [ p^ . 
The distribution of the eigenvalues and the eigenvectors, 
which are independent for the GOE or GUE, start to cou- 
ple when the disorder is increased and eventually break 
the O(A^) or U(A^) invariance. This breaking of the sym- 
metry appears in the large s behaviour of P{s) explaining 
why one can have a symmetry at the critical point which 
is independent of the magnetic field whereas the small 
s behaviour of P{s) is not. It has to be noted that the 
breaking of the O(A^) or V{N) invariance is at the origin 
of some interesting new results related with the lifting of 
the sum rule prohibition that is observed when studying 
the two- level correlation function jl^ . 

In summary, by means of the ELS method, we have 
studied the statistical properties of the energy spectrum 
of the Anderson model in a strong magnetic field at the 
MIT. For this we have considered the level spacing distri- 
bution P{s) as well as the Dyson- Mehta statistics, A3(L). 
Except for the small s behaviour of P(s), which depends 
on the presence or absence of a magnetic field but not 
on its strength, the other quantities studied, the large s 
behaviour of P{s) and A3(i) turn out to be independent 
of its presence or absence and thus of the time reversal 
symmetry. This suggests that the MIT, which is defined 
by the symmetry of the system at the critical point, is 



defined by a new universality (CE) class which is inde- 
pendent of the presence or absence of time reversal sym- 
metry. 
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